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1. INTRODUCTION 

The problem of electromagnetic wave generation by an electrodynamic tethered 
satellite system is important both for the ordinary operation of such systems and for 
their possible application as orbiting transmitters. The tether’s ionospheric “circuit 
closure" problem is closely linked with the propagation of charge-carrying 
electromagnetic wave packets away from the tethered system (Estes, 19881. 

Previous analyses of the waves generated by large conductors moving through a 
magnetoplasma (in our case a tethered system moving through the ionosphere) [Drell, et 
al, 1965; Belcastro, etai, 1982; Rasmussen, etal, 1985; Dobrowolny and Veltri, 1986; 
Barnett and Olbert, 1986; Estes, 1988] have considered the conductor to be Immersed in 
an infinite plasma medium. When (he boundary with the atmosphere is far enough 
away, this serves as a useful approximation for calculating the ionospheric waves and 
estimating their contribution to the tethered system’s electrical impedance; but it tells 
us nothing about the electromagnetic signal one should expect to be associated with the 
tethered system in the atmosphere or on the Earth’s surface. Heretofore there has not 
been a systematic treatment of the wave reflections and other effects of 
nonuniformities in the plasma medium. The inclusion of ion-neutral collisions in- 
troduces the possibility of coupling to other ionospheric wave modes. The work 
reported here represents a step towards a solution to the problem that takes into 
account the effects of boundaries and of vertical variations in plasma density, collision 
frequencies, and ion species. 

Section 2 recapitulates the theory of Alfven wave packet generation by an 
electrodynamic tethered system in an infinite plasma medium. A brief summary of 



Page 2 


previous work on the problem, including an assessment of where things stand with 
regard to controversial points, is given there. 

In Section 3 we present the generalization of our previous analysis that serves as 
the starting point — the incident wave-packet — for our treatment of the problem with 
boundaries and a nonuniform plasma. 

We present the physics on which our method of attack is based in Section 4. 
There we have examined the consequences of the presence of the boundaries and the 
vertical nonuniformity of the ionosphere. One of the most significant new features to 
emerge when ion-neutral collisions are taken into account is the coupling of the Alfven 
waves to the fast magnetosonic wave. This latter wave is important, as It may be 
confined by vertical variations in the Alfven speed to a sort of leaky ionospheric wave 
guide (centered at plasma density maximum in the F-layer), the resonances of which 
could be of great importance to the signal received on the Earth’s surface. We take the 
infinite medium solution for the case where the (uniform) geomagnetic field makes an 
arbitrary angle with the vertical as the incident wave-packet. Even without a full so- 
lution, a number of conclusions can be drawn, the most important of which may be that 
the electromagnetic field associated with the operation of a steady -current tethered 
system will probably be too weak to detect on the Earth's surface, even for large tether 
currents. This is due to the total reflection of the incident wave at the atmospheric 
boundary and the inability of a steady-current tethered system to excite the ionospheric 
wave-guide. 

An outline of our approach to the numerical problem is then given in Section 5. 
Given our conclusions about the weakness of the steady-current signal on the Earth’s 
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surface, we consider slowly varying tether currents.. We propose to use numerical 
integrations and boundary conditions consistent with a conducting Earth to obtain the 
solution for the horizontal electromagnetic field components at the boundary of the 
ionosphere with the atmospheric cavity. The proposed method, which is rather com- 
plicated, involves the use of Budden admittance matrices to take into account the 
incident and reflected Alfven waves and the fast magnetosonic wave to which the 
Alfven waves are coupled. We have begun software development for this project and 
plan to continue work on this problem under another NASA contract. 

A summary of our conclusions and planned future work is found in Section 6. 
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2. REVIEW OF PREVIOUS RESULTS ON TETHER WAVE GENERATION 

The theory developed in the prior Smithsonian Astrophysical Observatory 
(SAO) study of electromagnetic wave generation by an electrodynamic tether is 
reported in detail in the NASA report NAG8-551. A slightly modified version of a 
section of the report appeared in The Journal of Geophysical Research [Estes, 1988]. 
The publication of the results enables us to consider them in the light of the 
response they have received from other investigators in the field. 

Our study built upon previous ones [Barnett and Olbert, 1986; Belcastro et 
al., 1982; Dobrowolny and Veltri, 1986; Drell et al., 1965; Rasmussen et al., 
1985] but sought to obtain more believable results by using a more realistic model 
for the electric current distribution of the electrodynamic tether. Previous 
investigators had modeled the tether current system as a long cylinder. They had 
either considered an orbiting canister with a dimension along the direction of flight 
measured in tens of meters or, at the other extreme, an orbiting wire. Neither 
approach took into account the peculiar dumbbell shape of a tethered satellite 
system, which consists of a long, narrow wire terminated by satellites (or plasma 
clouds) with dimensions along the direction of motion that are much greater than 
the tether diameter. The SAO study showed that the orbiting canister approach, 
although it had not been thoroughly justified before, should provide a reasonable 
approximation for the purposes of calculating the wave impedance and ionospheric 
currents associated with steady-state (constant current) operation of an electrody- 
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namic tethered satellite system, provided that the current distribution assumes a 
constant value along the length of the tether, i.e., an insulated tether is assumed. 
This is true because it is primarily the changing fields associated with the injection 
of charge into the ionosphere at the ends of the moving system that drive the 
electromagnetic waves. Thus, any two moving systems that inject charge into the 
ionosphere at the same rate over the same area are equivalent from the standpoint 
of wave generation. By the same token, the orbiting wire model was found to be 
a poor model for an electrodynamic tether because it undervalued the dimensions 
of the terminating, charge-exchanging parts of the system by orders of magnitude. 

Barnett and Olbert [1986] of MIT had used calculations based on an 
orbiting wire model to advance the idea that the wave impedance associated with 
the operation of a constant current electrodynamic tethered satellite system would 
be on the order of 10,000 to 100,000 Ohms, as opposed to previous estimates on 
the order of an Ohm. This was not an academic question, since such high 
impedance values would have precluded the use of tethered systems for any 
applications that required substantial currents. This conflict has now been 
resolved, and Barnett and Olbert have acknowledged that their use of the orbiting 
wire model led to a great overestimation of the wave impedance associated with 
the frequency band lying between the lower hybrid frequency and the electron 
cyclotron frequency. In a private communication they expressed complete 
agreement with the SAO study on this point, and in a recent publication [Barnett 
et al., 1988] they credit R. Estes with pointing out that it is only the 
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system/ionosphere charge-exchange region that matters for wave excitation within 
the steady-state, cold plasma models considered up until now. 

Agreement has not yet been reached with M. Dobrowolny and his Italian 
colleagues on a number of points, which involve the functional dependence of the 
Alfven wave impedance on the Alfven and satellite velocities as well as the system 
dimensions. The SAO study, in agreement with all other studies but that of 
Dobrowolny and Veltri [1986] found that this impedance varies linearly with the 
Alfven speed, while being independent of the satellite velocity within the 
approximations used. Dobrowolny and Veltri, however, had found an inverse 
dependence on the Alfven speed and a quadratic dependence on the satellite speed. 
They have based their defense of these results largely on the thesis that the 
quadratic dependence on the satellite speed is much more reasonable physically, 
since the effect disappears (as it should) if the system is not moving. 

We will only sketch the argument here. The details of their argument 
and our rebuttal can be found in the Appendix. The Dobrowolny paper also 
contained a physical argument to explain the divergence of their results from those 
in the 1965 work of Drell, et al. We have already pointed out flaws in that 
argument [Estes, 1988]. The Italian investigators still have not addressed the 
question of why their results differed from all of the other prior studies as well, 
studies which obtained quite different results from those of Drell, et al. and which 
are in fact the relevant comparisons. Although Dobrowolny and his co-workers 
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have not emphasized this part of their results, the linear dependence of the Alfven 
wave impedance on the tether length that they found is as important a difference 
from the results of other investigators as the inverse dependence on the Alfven 
speed is. 

Basically, Dobrowolny and his colleagues maintain that the bulk of the 
plasma current is a “dc” current having nothing to do with wave phenomena. We 
argue that this amounts to saying that a current pulse traveling down a 
transmission line is a dc current. This approach appears fundamentally wrong 
when dealing with plasma currents, which are essentially wave phenomena. 

Our analysis has shown how the excitation of Alfven waves as the dominant 
form of radiation depends on the dimension along the line-of-flight of the charge- 
exchange region between the system and the ionosphere. It was the fact that this 
dimension is necessarily on the order of meters that enabled us to discard the 
higher frequency radiation that Barnett and Olbert were then arguing to be 
important. The vanishingly small dimension used by Dobrowolny and Veltri 
should lead to the same swamping of the Alfven band by the lower hybrid band 
radiation observed by Barnett and Olbert in their orbiting wire calculations, but 
since Dobrowolny and Veltri are only looking at Alfven waves they do not deal 
with this issue. 

Having failed to find an error in the straightforward SAO analysis, 
Dobrowolny and his colleague less advanced physical arguments about why the 


Page 8 


SAO analysis could not be right. They maintain that their result makes sense 
physically because there is an explicit quadratic dependence on the velocity of the 
satellite in their expression for the wave impedance, while the SAO result (they 
argue) does not go to zero as v x does. 

In the Appendix we demonstrate that our analysis also gives zero 
contribution to the wave impednace when the velocity is zero. The question is 
then: for a given non-zero satellite velocity, what is the proper dependence for the 
wave impedance on the satellite velocity? We have found that the wave 
impedance is very similar to that for a bifilar transmission line and have argued 
that this makes physical sense. If one considers a real bifilar transmission line, the 
property of the medium between the two wires (which are electrically insulated 
from each other along their length — as the magnetic field lines effectively are in 
the ionosphere) which enters into the calculation is the dielectric constant. Now in 
the low frequency approximation, which is applicable to the case of Alfven waves, 
the dielectric constant of the plasma (the diagonal, perpendicular components of 
the dielectric tensor that is) is in fact independent of the frequency. That is, it is 
independent of the satellite velocity, which is the only source of time variation in 
the problem under consideration. 

Thus the Alfven wave impedance of an electrodynamic tethered satellite 
system should be roughly independent of the satellite velocity. We can see no 
argument at all that it should be quadratic in v x , nor has anyone made a physical 
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argument for this particular functional form. We think this quadratic dependence 
found by Dobrowolny and Veltri is just the result of their having made an overly 
complicated analysis with some physically inconsistent assumptions. We demon- 
strate in the Appendix that, using the same starting point as the Italian 
investigators, we arrive again at our original results following two different routes. 

Having carefully re-examined our earlier results, we proceed to them as a 
starting point for the solution to the problem of the propagation of tether- 
generated ionospheric waves. 
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3. GENERALIZATION OF THE PREVIOUS RESULTS 

Modeling of the electrodynamic tethered satellite system has heretofore 
been based on some unrealistic assumptions regarding uniformity of the plasma 
medium and constancy of the tether current. In addition, the simplest relationship 
has been chosen for the system velocity vector, magnetic field vector, and the 
vertical. These three directions have typically been taken to define a three- 
dimensional orthogonal system. 

Taking into account the variations in plasma conductivity, density, collision 
frequencies, etc. in the ionosphere and the existence of the non-conducting 
atmosphere between the ionosphere and the earth will be the most difficult part of 
the process of obtaining a more realistic tether model. It will be an iterative 
project to solve the necessary boundary value problems. The analysis presented in 
this report is part of that project. Note that we use Gaussian units throughout 
our calculations. 

As a first step, we generalize our previous results on steady-current tethers 
to the case of arbitrary angles between the vertical direction, the satellite velocity 
vector, and the magnetic field vector. Let us consider the case for which the 
geomagnetic field is not perpendicular to the vertical, but the orbital velocity 
vector of the system is perpendicular to the plane of the vertical and the magnetic 
field (See Figure 3.1). This would correspond to a system at the maximum 
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excursion in latitude for a non-equatorial orbit. Having previously demonstrated 
the equivalence of an orbiting “ribbon” current distribution and the idealized 
dumbbell tether current distribution used in the earlier analysis (NAG8-551), we 
can conveniently define the tether current distribution as 

7 = W x ' - L x /2)~ H(x> + L x /2)] 

x 

. [H(yi - L/2)-H{yt +L/2)]6(z>) (1) 

where x 1 = x - v x t and y 1 lies along the vertical with y 1 = 0 at the middle of 
the tether. The y - axis, which is orthogonal to B (the z axis) and t; (the x axis) 
is indicated in Figure 3.1 as well. As in the previous analysis H(x) is the 
Heaviside function defined by 


H(x) =1, x > 0 

H[x) = 0, x < 0 (2) 


We now need k • Jk, where Jk is the Fourier transform of the tether current 
density. This is most conveniently calculated in the x,y f ,z> co-ordinate system, 
where 


k • Jk = kyljy^k 


( 3 ) 
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It is easy to obtain 


21 rf , \ sin k x L x /2 sink y ,L/2 

= v *("-H t— 


(i) 


Thus 


k • j = — 6 sm(k z L x /2 ) s j n ^^ ( ^/2) 

7T \ / 


( 5 ) 


Since 


7, = — k z (k- e) 
Z 4 7TW * V ' 


( 6 ) 


and 


47rio> ic • j 

C 2 (fc * 2 - W 2 £ x /c 2 ) 


( 7 ) 


have been obtained in the previous analysis, we can write the following expression 
for the Fourier transform of the Alfven wing or field-line sheet current: 
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Jz 


- 21 

it 


K 

(< K ? - w 2 e x /c 2 ) 
s\n(k z L x /2) 

k x L z 


. f( k v 

sin I - — 


cos 0 + k z sin 0)L 


) 


6(u - k x v z ) 


( 8 ) 


which differs from the earlier results for mutually orthogonal vectors only in the 
argument of one of the sine factors, which takes into account the dip angle of the 
magnetic field. 


The problem is now to obtain the inverse Fourier transform of (8). The 
integrals over u> and k y proceed very much as in the previous analysis, only with 
the results showing an obvious change due to the inclination of the magnetic field 
with respect to the horizontal. That is, the charge-exchange regions at the ends of 
the system are now located at y = ± (L/2) cos ft The integral over k z also shows 
the effect of the magnetic field’s inclination. The final integral over k x takes the 
form 


2nL, 

r Ko 


S^y - L cos 0/2^j 2H - L sin 0/2^J — 1 j 

J sm(Jc x L z /2) ex p^,^^ x f + fl/2)a] j dk x 

- —j- S^y + Lcos 0/2^ \2H[z + Lsin 0/2) - l] 

K 

J sto(k x L z /2) £ s i n 0/2)a]j dk z 


( 9 ) 
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1 

[i - (k x v x /n ci ) 2 ]H 

Although this expression looks more complicated than the corresponding expression 
from the previous analysis, the basic content is the same. The Alfven wings are 
seen to be field-line sheet currents at the ends of the system. The sign change in 
J z occurs at the charge-exchange interface as before, but these are now located at 
z = ± {L/ 2) sinfl. The top and bottom wings are connected by the condition of 
current continuity but otherwise they appear to be independent phenomena 
generated by the disturbances at their respective ends of the system. Except for 
the shift in lines of discontinuity in J z to coincide with those traced by the charge 
exchange terminals, the Alfven wing solutions are the same as before. The new 
solution clearly reduces to the old one when the angle 6 goes to zero. 

Another generalization of the previous analysis results when we allow for an 
arbitrary angle between the velocity vector of the system and the magnetic field 
lines. Let us consider the case when the tether lies along the vertical (y) axis and 
the magnetic field is in the horizontal plane antiparallel to the z-axis. The velocity 
vector, which lies in the horizontal plane, is allowed to be at an arbitrary angle <p 
with respect to the x-axis ((x,y,z) being a mutually orthogonal set of co-ordinate 
axes), as shown in Figure 3.2. 

As before, the dimension of interest is that of the charge-exchange interface 
perpendicular to the magnetic field lines between the system and the plasma. For 


where 


■■ ; 
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» 


Figure 3.2. Tethered system geometry (viewed from below) for velocity vector not 

perpendicular to magnetic field. 
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a spherical terminating satellite this would be the same dimension L x as before. 
The value of the current in the tether, assuming it to be induced strictly by the 
system’s motion, would be reduced by the factor cos <p since the induced voltage 
depends on the component of the velocity perpendicular to B. Our analysis 
is done in terms of the tether current I without calculating the value of I 
beforehand, but it is good to keep this difference in mind. 

Once again we make use of the equivalence between an orbiting ribbon of 
current and the system under study to write the tether current density as 

iy = J- - (osin v>)l) [AT(lf - 1/2) - + i/2)] 

• [H(H - 4/2) - H(x l + 4/2)] (10) 

where x ' = x - [v cos <p)t. 

This leads to the Fourier transform result 

k • J = — 6(w - k z v sin <p - k x v cos v?)sin(^^ sin(^ g /2) 
i r \ 2 / k z L x 


Our expression for the Fourier transform of J z is identical to the previous one 
((14) in our JGR paper) except for the changed argument of the delta function. 
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The integration over u> and k y can proceed as in the JGR article. The 
singularities on the k z real axis now occur at solutions of 

„\ 2 (k z cos<p+k z s\n<p) 

-) t 1 1 c (12) 

Va J (l - {k x cos <p + k 2 sin <p) 2 (v/n ci ) 2 ) 

From our previous results we can infer that, so long as <p is not too near tt/ 2, the 
condition 



k x cos <p | > | k z sin <p 


(13) 


should hold. 

Then we obtain 


K 


± 


k z v 

V A 



COS f 


k z v cos 



(14) 


for the singularities; i.e,, we just replace v z with ucos <p in our previous analysis. 
The other change occurs in the complex exponentials left in the integrand 


g xk x (x- vco$<pt) gik 2 (z - vsin<£>*) 
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After the integration over k z we obtain 


U X '<V’ Z ') = {[4(» - i/2)- 6(y + i/2)] 


[2ff(:')-l] ■ / 

Jo 


Kq 

dk 7 


sin(fc I 4/2) 


. e «M 2 ' +az ') | 


(15) 


where x ' = x - vt cosip, 

2 ' - - t;t sin£>, and 

_ V COSip 

V A\jl- (M COStp/n cl ^j 

Without getting all the details, we can see that this implies Alfven wings that look 
more or less like those shown in Figure 3.3. Seen from above it would appear as a 
winged structure flying along at an angle rather than “straight ahead.” 

There are a number of subtleties that would seem to be important, such as 
the actual shapes of the satellites, however, so that this analysis can only be 
viewed as a first approximation to an adequate description. 
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Figure 3.3. Alfven wing motion (viewed from below) for velocity vector not 

perpendicular to magnetic field. 
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4. ELECTRODYNAMIC TETHER IN A BOUNDED IONOSPHERE: 
BASIC PHYSICS AND METHODS 

Having generalized our earlier solution for the Alfven wings, we can now apply 
it to the case of a bounded, non-uniform ionosphere. In order to simplify the problem, 
we consider the tethered system to be immersed in a uniform geomagnetic field. The 
system’s motion is taken to be in a horizontal plane above the flat, conductive Earth. 
Its velocity vector is perpendicular to the plane formed by the vertical and the geomag- 
netic field lines. This corresponds to the situation described in the first part of the 
preceding section and illustrated in Figure 3.1. We will adhere to the notation of the 
previous section in the following analysis. 

The first significant new feature we introduce to the problem is the presence of 
boundaries: the one between the ionosphere and the atmosphere and the one between 
the atmosphere and the Earth. For simplicity, we begin our analysis with the atmo- 
spheric cavity considered as a vacuum and the Earth as a perfect conductor. For elec- 
tromagnetic problems the first approximation is reasonable, and the second approxi- 
mation greatly simplifies the boundary value problem at the Earth. If the ocean sur- 
face is considered, it is justifiable as a first approximation since its main consequence 
is a small horizontal electric field component at the surface. 

A complete, self-consistent solution to the problem of an electrodynamic teth- 
ered satellite system operating in the more realistic environment described above is 
beyond the reach of this analysis. We make the assumption that the tethered system is 
sufficiently far from the atmospheric boundary (or any steep gradients in plasma pa- 
rameters) that we need not be concerned with the boundaries’ effects on the system. 
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That is, we assume that the infinite-medium solution previously obtained is a reason- 
able approximation to the “incident” wave-packet generated by the system operating 
in the bounded ionosphere. We do not claim to have established criteria by which to 
judge what “far enough away" means, but our feeling is that the approach can be fruit- 
ful if the system is sufficiently far from a boundary that it (including its local Alfven 
wing structure) moves into a region previously unoccupied by any part of itself within 
twice the transit time of an Alfven wave to the boundary; i.e., if the system has “moved 
on” before the reflected wave arrives. It should be noted that TSS- 1 cannot be expected 
to satisfy this condition most of the time because of the steep plasma density gradients 
between its 300 km orbit and the atmospheric cavity. For now we assume a tethered 
system high in the ionosphere. We restrict ourselves to the steady-state operation of 
such a system and take the solution of the previous section as the incident Alfven wave 
packet. We allow only vertical variations in ionospheric quantities. 

The complexity of our problem, even in the simplified form stated above, re- 
quires a numerical analysis. There are. however, a number of observations that can be 
made based on the fundamental physics of the system under consideration. Our ap- 
proach follows the general outlines of the analyses of ionospheric waves made by P. 
Greifinger [1972), C. and P. Greifinger 11973], and Rudenko, et at, [1985]. The particu- 
larities of our moving source require some modifications to the analysis from the out- 
set. however. In order to utilize the formalism of the above-mentioned authors, we 
seek an Incident wave solution written in terms of plane waves in the horizontal 
plane, the amplitudes of these wave components being dependent on the vertical co-or- 
dinate, y' in our notation. A number of transformations, which will be defined in de- 
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tail below, are required to obtain this form for our incident wave-packet. We must also 
take into account the relationship that exists between the x component of the wave 
vector and the frequency (as seen in the plasma, i.e. the terrestrial, rest frame) for the 
steady-state operation of an electrodynamic tethered system: namely, the Doppler re- 
lation 0) = k x V x . This consequence of the constant tether current (which, it should be 
remembered, is an assumption about the nature of the interaction between the tethered 
system and the ionosphere) turns out to be of great importance to the electromagnetic 
field on the Earth’s surface, which we are hoping to calculate. At this initial stage of 
our analysis we restrict ourselves to the Alfven region k X V x < < , which is 

consistent with our consideration of systems in the upper regions of the ionosphere. 

Following the pattern used to derive expressions (4) and (8) in the preceding sec- 
tion, we obtain the following expressions* for the Fourier transformed electric field 
components 


E 


*.y 


SicoI_ 

k\c 2 


k x v x ) 



x,y 


(16) 


* In all that follows we will use the notation k y = k , etc. to represent vector 
components in the primed system of co-ordinates. 
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where the _L subscript refers to components perpendicular to the magnetic field, 
which lies along the negative z direction. 

Thus Maxwell's equations give us 




(17) 


Since E z — 0 by assumption. Maxwell’s equations further yield 


„ k z c „ 

B = + — —E 

y ,x - CO x , y 


(18) 


These equations represent only the Fourier amplitudes of the field compo- 
nents. The full expressions include the plane wave complex exponential factors 
exp (- iicot - k x)) . 

Our first step in obtaining the desired expressions for the horizontal field com- 
ponents in terms of horizontal plane waves is to carry out the inverse Fourier trans- 
form integration over k z . Since we are considering downward moving waves, which 
correspond to the negative z direction in our co-ordinate system, we close the contour 
of integration in the lower half of the complex k z plane and pick up the contribution 


of the pole at 
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k z 



( 19 ) 


We thus arrive at 



where k is evaluated with k . 


CO 


We now use 


k y y + k z z = k’ y y'+ k' z z' 


( 21 ) 


and the transformation of variables 


1 1 \dcodk x dk y 


= sin 


0 1 1 \dcodk x dk' z 


( 22 ) 
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to obtain 


E x = 


A | J jdcodk x dk' z {exp(- i[(Qt - k x x - k y y' - k Z z'~\ ) 


nL x c> 


sin 


S( co - k x v x ) 


n sin 


, 2 


( 23 ) 


E y = 


— A 2 — J J jd(odk x dk' z {cx p(- i\_cot - k x x - k' y y' - k' z z'~\ ) 

t±, x c 


sin 


S(a> - k x v x )- 


k x L x 


k\k x 




y cos 0 - k z sin 


( 24 ) 


where 


k 


y 


co 


. - - k 2 cot 0 
v , sin 0 2 


( 25 ) 


and 


k x = k x + (k y cos 0 - k z sin 0 ^ 


( 26 ) 
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Figure 4. 1 . The geometry of the problem 
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All vectors lie in the horizontal plane. S is parallel 
to the horizontal k vector. f3 is chosen to make a right- 
handed orthogonal system of (S 7 y', p) 


Figure 4.2. The co-ordinate system defined by the horizontal wave 
vector 
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Equations (18) and (19) imply 


B 


x,y 



(27) 


The horizontal components of the field are now at hand. The x-component is 
one horizontal component. The other horizontal component of the wave field (the z' - 
component) is easily obtained utilizing E z = B z = 0 


E? z = - E y sin 9 


(28) 


The corresponding expression holds true for the B field components. 

Following the approach of the references mentioned earlier (particularly 
Rudenko, et at), we now need to obtain the horizontal wave field components parallel 
and perpendicular to the horizontal wave vector. Figures 4. 1 and 4.2 show the co-ordi- 
nates to which we refer. Note that each horizontal wave vector (k x , k z ) defines a 
different co-ordinate system (8, fi) . This implies that we will have to carry out our 
numerical integration separately for each (k x , k z ) pair. In all that follows the 
vector k with no subscript refers to the horizontal wave vector of magnitude 
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k = 


J 


k\ + *7 


(29) 


In the 


( S, p) system the wave electric field components are 


E k 




k 


(30) 


E 


p 




(31) 


The corresponding expressions hold for the components of the B field. 

Having obtained the components of the Incident wave packet , we can consider 
some of the general physical characteristics of the system. The equations to be inte- 
grated in the ionosphere will be presented in detail in the next section. We will allow 
for vertical variations in plasma quantities. These variations include the sharp 
changes at the boundaries in our problem. Since the only variations in our model are 
in the vertical direction, the dependence of our incident Alfven wing wave packet com- 
ponents on the horizontal co-ordintates remains the same throughout and is given by 
the horizontal plane wave factor. Thus each component of the incident wave packet 
will arrive at the boundary between the ionosphere and the atmosphere with the same 
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horizontal wave vector that it had high in the ionosphere. This has important 
consequences for the solutions in the atmosphere and on the Earth. 

The equations for the horizontal field components in the atmospheric cavity 
are 


d *E 


5. P 


dy' 


+ 


(V . 2 

2 k 

V C 


E s,p=° 


(32) 
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5, P 


dy' 


r K-k 2 )B, =0 


V c 


S, P 


(33) 


where k , the magnitude of the horizontal wave vector, has the same value as in the 
ionosphere. Similarly, the frequency is unchanged across the boundary. 

The incident wave components contain the factor S( (0 - k x v x ) , a conse- 
quence of the steady-state operation that we have assumed up until now. An immedi- 

(-- k 2 ) 

ate consequence of this is that the factor y c 2 J in the second term of equations 
(39) and (39) is always negative, since \k x v x | < | kc\ always. 

This means that the Alfven wings generated by the steady state operation of an 
electrodynamic tethered satellite system will not propagate into the atmospheric cav- 
ity. That is, there is total reflection at the atmosphere's boundary with the iono- 
sphere. Our solution corresponds to a surface wave at the ionospheric boundary. 
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Taking into account the perfect conductor boundary condition at the Earth’s 
surface, we obtain the solution in the atmospheric cavity 


B s. p p( y b ) 


sinh(p(/-(/ t -g))) 
si nh(p H ) 


(34) 


B r 


jrcE e (y>) 


cosh (p(/-( y\ - H))) 
si nh(pH ) 


(35) 


B S=~ 


l p c 


CO M4) 


cosh (p( /- ( y\ - H))) 
sinh(p H ) 


(36) 


where 



and y' b and H are the values of y' at the ionospheric boundary and the distance 

of this boundary from the Earth’s surface , respectively. 

Equations (35)-(36) show that the ratio of the magnetic field on the Earth’s sur- 
face to that at the ionospheric boundaiy is 
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^s. p( y b) 


= (cosh( pH)) 


(38) 


The height of the ionosphere may be taken to be around 100 km. The conse- 
quence of equations (35)-(38) is that the image of the Alfven wings on the Earth’s sur- 
face will be much wider than the wings are in the ionosphere, since only the wave- 
packet components with horizontal wavelengths of hundreds of kilometers will escape 
severe attenuation. Since such long wavelength components make up only a small 
fraction of the wavepacket for a reasonably sized tethered system and since the noise 
level is much higher for the lower regions of the ULF band, we tentatively conclude 
that the magnetic field image of the Alfven wings on the Earth’s surface will probably 
be too weak to detect in the case of a steady-current tether, even one with a high cur- 
rent. 

Hughes and Southwood [1976] reached similar conclusions about the 
“shielding’’ of ionospheric disturbances with short horizontal wavelengths. These au- 
thors emphasized the role of Hall currents in the lower part of the ionosphere in 
reducing the B ^ component of the ionospheric waves; but it is obvious from the 

analysis that the result is quite general since equation (35) is a consequence of 
Maxwell’s equations and the boundary conditions and would hold for different models 
of the ionospheric conductivity 

So far we have not discussed the possible effects of vertical variations in 
plasma density and ionic composition. Since the Alfven speed depends directly on 
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these quantities, their variation changes the effective dielectric constant of the 
plasma. This has consequences for the waves that can propagate in the ionosphere. 
The dispersion relation for the anisotropic Alfven wave is 


co= k 2 v A 


(39) 


The dispersion relation for the fast magnetosonic (isotropic Alfven ) wave is 


2 , 2 2 
0) =k tol V A 


(40) 


where k u is the total wave vector, including horizontal and vertical components. 
For a given horizontal wave-vector, frequency, and Alfven speed combination it may 
be impossible to find a real vertical wave vector component that satisfies the disper- 
sion relation(40). In this case the fast magnetosonic mode is evanescent. This is the 
case for the steady-current electrodynamic tether, which excites waves satisfying (39), 
but by virtue of the CO = k x v x condition cannot satisfiy (40) for any real vertical 
wave vector component. The fast magnetosonic mode has thus been discarded in our 
calculations of the Alfven wings. As an aside, we note that this neglect of the fast mag- 
netosonic mode, quite justifiable in the infinite plasma case, where only propagating 
waves are considered, should probably be re-examined in the bounded plasma case. 

We have seen that for a given horizontal k vector and frequency a, the isotropic 
wave will propagate or not depending on the Alfven speed. The Alfven speed has a 
minimum value in the F-region of the ionosphere, so an isotropic wave can be con- 
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fined to the regionaround this minimum — the ionospheric wave guide. Figure 4.3 
schematically depicts the ionospheric wave-guide. The Greifinger and Rudenko refer- 
ences discuss the ionospheric wave-guide and note the occurrence of wave-guide 
resonances. Rudenko, et al. make the point that an Alfven wave, incident from high 
in the ionosphere, can couple to the fast magnetosonic wave , through ion-neutral 
collisions, within the ionospheric wave-guide region for the particular frequency and 
wave-vector combination. 

The ionospheric wave-guide traps electromagnetic energy in a horizontal layer 
of the ionosphere. Stimulating resonances of the ionospheric wave-guide with an 
electrodynamic tether would seem to be a possible way of overcoming the difficulties 
in obtaining a measurable signal on the Earth's surface. It is easy to see that no 
steady-current electrodynamic tethered system can achieve this in the Earth’s iono- 
sphere, however. Since an orbiting steady-current tether cannot stimulate propagat- 
ing fast magnetosonic waves even in the region of minimum Alfven speed, it is impos- 
sible for Alfven waves generated by a steady-current tether anywhere in the iono- 
sphere to excite propagating fast magnetosonic waves. There is no ionospheric wave- 
guide for steady -current tethers. 

The next phase of our research into the problem of tether-generated 
electromagnetic waves will focus on the ionospheric wave-guide and its possible 
excitation by electrodynamic tethered satellite systems with time-varying currents. 
In the next section we map out the method we plan to use in our numerical 
investigation of this question. 
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Figure 4.3. The ionospheric wave-guide and other layers of the medium 
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5. METHOD FOR INVESTIGATING WAVE-GUIDE EXCITATION 

We should state from the outset that our approach is not applicable to all types 
of tether wave problems. It should be possible to extend the method, however, once it 
has proven its value. In its present form, the method can be applied to the case where 
the tethered satellite system is well above the ionospheric wave-guide for the the waves 
it generates. That is, it applies to a system that directly excites only anisotropic 
(shear) Alfven waves. The integration method would thus apply to the case of a steady- 
current tether; but, as we have already noted, the ionospheric wave-guide concept does 
not apply to such systems, since there is no region of the ionosphere that will propa- 
gate isotropic Alfven waves that concurrently satisfy the shear Alfven wave dispersion 
relation appropriate to a steady-current tether. Thus we would expect the method to 
have its greatest utility in the case of slowly-varying tether currents with frequency 
components such that coupling to ionospheric wave-guide modes occurs. 

It should also be noted that we are talking only about variations in the tether 
current sufficiently slow that all of the previous work on steady-current tether wave 
generation applies with only slight modification. This means that the tether current 
distribution can still be considered as independent of the vertical co-ordinate along 
the tether’s length. The tether is not functioning as an antenna in the normal sense of 
the word in this case, any more than it does in the steady-current case. It is still a 
source of varying current to the ionospheric transmission line, a concept that we have 
emphasized in our previous analysis. The ionospheric waves generated by the tether 
are due to the charge it injects into the ionosphere at each end of the system, rather 


than to direct radiation from the tether. 
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We assume a tether current distribution that is just expression (1) multiplied by 
cos (cojt ) , where G) 0 is the oscillation frequency of the tether current. How this oscil- 
lation is achieved need not concern us too much at this stage. Let us suppose that the 
tethered system is equipped with a suitable power supply to enable it to reverse the cur- 
rent. Any time variation can be broken into its Fourier components, so the choice is 
not restrictive, except that we are explicitly excluding spatial variations in the tether 
current along its length, which implicitly limits the range of frequencies we can con- 
sider. 

Expressions (23) and (24), which define the incident wave packet, are modified 
only by the replacement of the delta function 5 (g) — k x V x ) by 

( 5(g> + o) 0 - k x v x )+ 5(g)- c o 0 - k x v x )) 
2 < 41 > 

Having derived the incident Alfven wave-packet, we now turn our attention to 
the numerical methods we propose to use for solving the problem. We are following 
the method outlined in Rudenko, et aL, However, our notation differs in some respects 
from that used in the Rudenko reference. In addition there are a few serious typo- 
graphical errors in that paper. For these reasons, and to make this report stand alone, 
we will outline the method here. 

The incident wave solution we have obtained applies to a region where the ion- 
neutral collision frequency is negligible and the Alfven speed is constant. Anticipat- 
ing the variations in plasma quantities that occur lower in the ionosphere, we express 
Maxwell’s equations in the following compact way : 
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_ .cLF 
1 dy' 


= G F 


(42) 


where we define F as the four component vector 


F=($,E IJ H s ,H p ) 


(43) 


with the third and fourth components given by 


H s 



(44) 



(45) 


For our numerical work we use dimensionless quantities. The dimensionless 
frequency is defined by 



( 46 ) 
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where X is a linear scaling factor and v Q is the Alfven speed in the region in which 
the wave-packets originate. The variable y' is just the vertical spatial co-ordinate in 

units of X . 

In expression (42) G represents the 4 X4 matrix 


G = 


Si S 2 ^ 
S 3 J 


(47) 


defined by the 2 X2 submatrices 


Si 


U 2 cotd - k x cotd N 
0 0 , 


(48) 


•■(if) 


(49) 


' 1 ( r 2 r 2 2 N 

U - U x cos 2 0 


g 3 = 


v 


V 


ft sirfO 


— 2 ( 
CO 
_ 2 
V 


_2 t 

CO 

- 2 
V 


- 2 
0 1^2 


}i x k x Ti z cot 6 ji 2 


V 


k 


sinB 


V 

— 2 
CO 

— Ml 

V 


f 2 2 

£ - % v cos 2 8 


sinO 

\ 




. 2 

k sin 0 


- I 


(50) 



Page 41 


r - k z cotd 0 > 
k - k x cotO 0 j 


In the expressions above k z and k x refer to the dimensionless horizontal 
wave-vector components defined by the scaling factor X ■ The normalized Alfven 

Za 

speed V is defined by u o . 

The quantities 


^1 = 



( 52 ) 


and 


0 )]^ ( 53 ) 

along with the Alfven speed, express the dielectric properties of the ionospheric 
medium. Here v { and Q ci are the ion-neutral collision frequency and the ion cy- 
clotron frequency, respectively. The Alfven speed . V , and Q ci all vary with altitude 


in our model. 
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For the frequency range and horizontal wave-vectors to which we are limiting 
our analysis, there are four well-defined independent solutions to equation (42) in the 
upper ionosphere: two shear Alfven solutions, corresponding to positive and negative 
vertical wave-vector components, and two isotropic Alfven solutions, one that grows 
exponentially with increasing y' and another that falls off exponentially with in- 
creasing y‘. The shear Alfven solution with negative vertical wave-vector component 
corresponds to our incident wave. The amplitude for each horizontal wave-vector 
component is given by (23) -(27), as modified by the substitution of (41) as specified 
above. The upward traveling shear Alfven solution then corresponds to a reflected 
wave. Of the two isotropic Alfven wave solutions, only the upwardly decreasing solu- 
tion makes physical sense. This corresponds to leakage of the ducted fast magne- 
tosonic wave from the ionospheric wave-guide. Thus we are left with three physically 
meaningful solutions to equation (42) at “infinity” (which is the short-hand expres- 
sion we will use to indicate the location of our tethered system, high in the iono- 
sphere). 

While we know the amplitude of the incident wave-packet solution at infinity, 
the amplitudes of the other two solutions are unknown and must be the result of com- 
plex interactions in the lower ionosphere. It turns out that knowledge of the func- 
tional form of these solutions at infinity, combined with the boundary conditions at 
the ionosphere /atmosphere interface and our knowledge of the functional form of the 
solutions in the atmospheric cavity (Equations (34)-(37)) will suffice to determine a so- 
lution at the boundary between the atmosphere and ionosphere and, hence, on the 


Earth’s surface. 
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Let us now sketch the means by which this can be accomplished. First we in- 
troduce the admittance matrices, which are variations of those defined by Budden in 
the “Bible" of ionospheric wave physics The Propagation of Radio Waves [1985]. We 
define the admittance matrix A (co, k, y') as the 2 X2 matrix that satisfies the 
equation 


( 


V 


oF\ + 0F* ' 

< + PFi j 


= E, y') 

V 


aF\ + fi F J X ' 

<*2 + j 


(59) 


where F ‘ and are two solutions of (42) and a and P are two arbitrary complex 
constants. In other words, the admittance matrix transforms a linear combination of 
the electric field components of the two independent solutions into the same linear 
combination of the corresponding magnetic field components, as defined by equations 
(43M45). 

This is equivalent to 


( 


T J ' 
3 

< j 


A ( y') 


F\ F 
F* F. 


J \ 


2 J 


(60) 
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It is straightforward to obtain from (42) and (60) the following differential 
equation for an admittance matrix A defined as in (59) and (60): 


- l 


dA 

dy ' 


= - AA- A g,+ g 4 A +g. 


(61) 


where the£/ are the matrices defined in (48)-(51). 

The admittance matrices contain the ratios of electromagnetic field compo- 
nents rather than their absolute values. This has an important advantage for numeri- 
cal integrations down through the ionosphere, since it avoids the problem of numeri- 
cal swamping brought on by the exponential growth of an initially small mix of the 
solution that grows as the altitude decreases. 

Since we know the functional form of the solutions at infinity, we can con- 
struct the admittance matrices at infinity in the following way: 


d.j) 
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*3 *3 

xi x J 
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x* x J 


j ^ 


-1 


2 J 


(62) 
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where the X* are the known functional forms of the F * solutions at infinity. 
These can then be taken as the initial values for numerical integration down to the 
boundary with the atmosphere. 

A the boundary matrix at the atmosphere/ionosphere boundary is defined by 


(H n)i 

K( r.)J 


= A 
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(H y\) 

w n) 


\ 

J 


( 63 ) 


where the F corresponds to the total solution at the boundary, including the contribu- 
tions from the reflected and ducted waves, as well as the incident wave. 

For the case of a perfectly conducting Earth this corresponds to 
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0 , 

= iAcoth ( pH ) 

\ c P) 
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( 64 ) 


The solution at the boundary may be written as 


F a? t ) = B,F '( y b ) + B 2 f\ y\) + B,F\ y\) 


( 65 ) 
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1 2 

where F and F denote the reflected and incident shear Alfven solutions, respec- 
tively, and F denotes the ducted wave solution, with the B i complex constants. The 
coefficient B 2 is known from the incident shear Alfven solution at infinity. 

An obvious consequence of the definition of the admittance matrices (59) is 

that 


/) - A (k ' J) ( to, £, /)) 




= 0 


( 66 ) 


for i * k . 

We can utilize expressions (59), (63), (65), and (66) to obtain the following equa- 
tion, which is true on the boundary of the ionosphere and the atmosphere: 


(a 6 (o), E) - A ( 1 ’ 3 ) (o),^,/J)^ y' ^ 

( A (2 ' 3) ( w, h, y'J - A (U) ((0, k, y' b )) B 2 


(67) 


f (2)\ 

F* 

( 2 ) 

V F 2 J 


This expression, which relates the electric field components of the total solu- 
tion (65) on the boundary to the electric field components of the incident wave solution 
on the boundary by means of the admittance and boundary matrices, is the basis for 
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our numerical method. We can invert the matrix ( A (to, *)- A (1,3 ) (w, E,/ 6 )) 
on the left hand side of (67) to obtain the electric field components on the boundary. 
Then the boundary matrix A& yields the magnetic field components on the boundary. 
The atmospheric cavity solution (35)-(36) yields the magnetic field on the Earth’s 
surface and anywhere in the atmosphere. The required admittance matrices and the 
incident wave solution on the boundary are to be obtained by means of numerical 
integration of the equations (61) and (42). 

All of the analysis outlined in the preceding paragraphs must be carried out for 
each horizontal wave vector component. The complete solution on the Earth’s surface 
is obtained by summing over all these solutions to obtain the inverse Fourier trans- 
form. A sort of flow chart of the method is shown in Figure 5.1. 

This is the scheme we propose to carry out to obtain the first approximation to 
the electromagnetic field on the Earth’s surface due to an orbiting electrodynamic 
tethered satellite system with a slowly varying current. It requires a model for the 
ionosphere that includes vertical profiles of the plasma density, the effective ion 
mass, and the ion-neutral collision density. Other assumptions of the theory are dis- 
cussed in the concluding chapter. We have begun software development to carry out the 
program of numerical analysis. 
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Figure 5.1. Flow chart of the complete 
computation 
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6. CONCLUSIONS AND PERSPECTIVE 

During the course of this work we have re-examined our previous analysis of 
ionospheric wave excitation. by electrodynamic tethers in an infinite medium. Our 
rebuttal to the criticisms made by Dobrowolny and his co-workers is contained in the 
Appendix. We remain convinced that our analysis was sound, and it appears that most 
investigators in the field concur. We note that Barnett and Olbert of M.I.T. have ac- 
knowledged that our analysis, as contained in the JGR paper, was correct in pointing 
out that their orbiting wire model of the tether lacked physical sense and led to a gross 
overestimation of the contribution of the lower hybrid frequency band to the tether 
wave impedance. 

We have generalized the analysis to the case where the geomagnetic field lines 
make an arbitrary angle with the vertical. The generalized solution thus obtained has 
then served as the starting point for a first approximation to the much more difficult 
problem of tether-generated Alfven wave packet propagation through a non-uniform 
ionosphere bounded by the atmosphere, which is in turn bounded below by a conduct- 
ing Earth. 

We have discussed the physics of the ionospheric wave-guide for fast magne- 
tosonic waves and noted the possible importance of the excitation of ionospheric 
wave-guide modes coupled to the incident shear Alfven waves generated by the tethered 
system. A steady-current electrodynamic tether will not excite the ionospheric wave- 
guide. however, as we demonstrate in Section 4. 

We have shown that it is possible to gain considerable insight into the physics 
of an electrodynamic tether in a bounded ionosphere without solving the complete 
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problem. In addition to the result stated in the previous paragraph, it was also possi- 
ble to show that the Alfven wings generated by a steady-current tether are completely 
reflected at the boundary of the ionosphere with the atmosphere. Thus there is no di- 
rect propagation of the waves into the atmosphere, only a leakage of the associated 
field. On the Earth’s surface only wave components with horizontal wave-lengths in 
the 100 km and greater range would have the same strength as at the ionospheric 
boundary . Thus the image of the Alfven wings on the Earth’s surface would be much 
wider than in the ionosphere. The signal strength would probably be too weak to de- 
tect, even for high current values. 

Thus the problem of time-varying tether currents assumes even greater impor- 
tance from the standpoint of the signal reaching the Earth. We have made the initial 
steps in this analysis, carrying the generalization one step further to consider the case 
of slowly varying tether currents. We have presented in some detail a scheme for ob- 
taining the magnetic field on the Earth’s surface for an electrodynamic tethered satel- 
lite system with a varying current in orbit above the peak in ionospheric electron den- 
sity, i.e., above the center of the ionospheric wave-guide. This method, which involves 
the numerical integration of the initial Alfven wave-packet and Budden admittance 
matrices down through the ionosphere, will be applied in our ongoing work on this 
problem. 

It is always good to keep in mind the limitations and assumptions of a theory. 
The ionospheric plasma is considered to be a linear, dielectric medium. As in previous 
analyses, we have assumed that the plasma temperature is not important to the waves 
under consideration, though we have included ion-neutral collisions in our new model 
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for the dielectric tensor. Based on previous experience, we have simplified the prob- 
lem by considering only wave components satisfying the magnetohydrodynamic cri- 
terion. Vertical variations in plasma density, ion content, and collision frequency are 
allowed: but horizontal variations, which could be important in practice, are not in- 
cluded. The Earth is flat, the magnetic field uniform. Furthermore, the incident 
wave-packet has been taken from the infinite medium solution. All of these assump- 
tions could eliminate effects that will turn out to be important. The extension of the 
analysis to higher frequency tether current variations is clearly a necessary step, but 
one that will require greater care than the analyses we have seen so far, since deter- 
mining the current distribution within the tether becomes crucial and requires a self- 
consistent solution to the problem of the interaction between the tether and the iono- 
sphere. 

While acknowledging that much work remains to be done, we feel that the work 
reported here represents an advance in our understanding of the problem of electro- 
magnetic wave generation and propagation from a tethered satellite system and points 


the way to further progress. 
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APPENDIX: An Examination of the Dobrowolny/Veltri Results and 

Criticism 

In this appendix we will consider in detail the analysis of Dobrowolny and 
Veltri, since Dobrowolny and his co-workers have continued to maintain that it 
alone has found the correct solution to the problem of tether wave generation. 

We want to make clear that as regards their criticism of our results we will 
be summarizing verbal remarks Dobrowolny and his colleague less have made 
rather than referring to anything they have written. We believe we have 
understood their criticisms and will attempt to make a fair presentation of them, 
while stating from the outset that their arguments have not convinced us. We 
have in fact rederived our results using their formulation of the problem. 

First let us consider their original analysis (as presented in the Nuovo 
Cimento article by Dobrowolny and Veltri) and compare it with the SAO analysis. 
Dobrowolny and Veltri use a long orbiting cylinder with square cross-section to 
model the tethered system. Since D, the tether’s transverse dimension, is stated to 
be on the order of 1 mm, it is clear that this is just another variation of the 
orbiting wire model, the inadequacy of which we have previously demonstrated. 
However, the authors do not deal at all with the issue of radiation in other 
frequency bands raised by Barrett and Olbert. Following a different approach 
from that of the SAO study, they begin by writing down the most general, all- 
inclusive formula for the power radiated by a current distribution in a plasma 
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(equations (9)-(12) in the original): 

P = T ~ j I % K ( 4 . 1 ) 

where 

A* = Tr\ { j (-A.2) 

and A,y is the co-factor of the plasma tensor Ay defined by 

Ay = n 2 ^ + Sjj'j + eij(k,u>) (A.3) 


using standard notation. Finally A — det Ay, and ey is the polarization vector 
defined by 


c. = 



(4.4) 


in terms of an “arbitrary complex vector a.” J( e ) is the external current, which 
we have denoted by j. Neither at this point nor at any later point do the authors 
ever specify the form of £,y. They maintain a high level of abstraction. The 
meaning of the polarization vector seems particularly obscure to us in this context. 
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For formula (A.l) the authors refer to Plasma Astrophysics: Nonthermal 

Processes in Diffuse Magnetized Plasmas, Volume 1 by D.B. Melrose. In fact the 
expression does not appear precisely this way in the reference. What Melrose 
writes is (equation (3.8) p. 65) 

' - f 7 / £ / <0 (A*) 

where we have slightly altered the notation to be consistent with Dobrowolny et al. 
The integrand in Melrose is the dot product of the Fourier transform of the wave 
electric field vector and the complex conjugate of the Fourier transform of the 
external current density, which is the tether current in our case. 

Before continuing with our discussion of the analysis of Dobrowolny, et al. 
let us apply this equation, i.e. the equation as it actually appears in Melrose, to the 
SAO results. 

We have already seen that the contribution of the j z component to the 
radiated power is small, so we consider the contribution of j y *E y . Referring to the 
previous SAO results, we obtain 
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Using equation (2.12) (p. 27) of Melrose we can make the substitution 


[6(u; - A^t/*)] 


2 


lim 

T — > oo 



Mz) 


(A.7) 


Then we obtain the expression for the power going into waves 


4/ 2 t f d(jjd 3 k u ~ k x v x )s\n 2 {k y L/2)sm{k x L x /2) 
Z z 7T 2 C 2 J kj 2 • (fc z 2 - W 2 6 x /c 2 ) 


(A.8) 


where factors of 27T have disappeared because we use a different convention from 
Melrose’s in defining the Fourier transform. Performing the integration over u 
and k z exactly as before, we get 


- ^ r «. i\ /, _ ( ,. 9) 

tvL z c 2 J_ Ko k ± 2 k x y \^c*y 


Utilizing the formula 
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/ sin 2 ax 

x 2 +j3 2 


dx 


Tp 

fa > 0, fie/S > o) 


M-10) 


found in Gradshteyn and Ryzhik’s tables (p. 447), we finally obtain 



(All) 


Dividing by I 2 to obtain the impedance Z^, we find we have arrived at an 
expression identical to the one derived previously in the SAO study. This exercise 
has demonstrated that by starting from the same point as Dobrowolny et al. we 
still arrive at the same result as previously. Since expression (A.5) is completely 
general, this had to be, assuming we had made no errors. 


The purpose of this exercise has been twofold: first to point out that if 

there are errors in the SAO analysis they must have occurred in the derivation 
of the expression for Since this derivation is easily followed, it should be 

possible to discover any errors made there. The point that follows from this is 
that we should be doubly careful in evaluating the analysis of Dobrowolny et al. 
which supposedly uses the same starting point. 


The SAO analysis made use of the E z = 0 approximation to good purpose. 
When this is not done, a great deal of unnecessary apparent complexity remains, 
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with a concurrent increase in the probability that errors in calculations or physical 
reasoning will occur. 

Let us turn our attention to equation (A.l), the starting point for the 
analysis of Dobrowolny and Veltri. We have already mentioned that it is itself a 
derived expression not found explicitly in the Melrose reference. Although the 
authors have not demonstrated the intermediate steps, we can attempt to 
reconstruct their chain of reasoning. 

Melrose gives the general equation ((3.16), p. 67) 


Ei{k, w) 



A(k,u) 



{A.12) 


This follows from the definitions and the basic wave equation. 
Inserting this into equation (A. 5) gives 


- lim 4ni f du f (fik Jj_ ^ 

T -* oo T J 2ttu> J (2tt) 3 A 


(A.13) 


Now for a given mode a, Melrose defines the polarization vector in the way 
Dobrowolny and Veltri seem to be doing for all modes, inclusively, in expression 
(A.4). We have not examined this issue thoroughly, but it is not immediately 
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obvious to us that this extension makes sense, i.e. that a single arbitrary complex 
vector a would yield physically meaningful results for all modes. Granting this 
point (They never say what a they use.), we then obtain 


A 




Ass 


(A.14) 


which leads to 


- lim % 

T -» oo 47r 3 T 




(A.15) 


This becomes identical with equation (A.l) if we make the substitution 


1 

d(fc, w) 


2 7T* 5 



(A.16) 


At this point, the meaning of 5(d) is somewhat ambiguous. Although the authors 
present the next step in their results without any explanation beyond a 
characterization of the intermediate calculations as “lengthy”, we gather from the 

factor | dX/du | 1 appearing in their result that they are taking the 5(d) to mean 

that A is a function of u, which depends parametrically on k insofar as the 
integration over u is concerned. Thus the integration gives a sum of terms from 
the poles of d _1 (i.e. the zeros of A). These poles would correspond to the 
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solutions of the dispersion relation A(ut,k) = 0. Therefore 



- u„(kfj/ 


dA_ 

dui 


W=U*t(*) 


where the sum is over the different wave modes. 

Since they were only going to consider Alfven waves, the authors could have 
saved themselves a lot of computational labor by continuing to follow the analysis 
presented in Melrose’s book. On page 56 of that work we find that following 
results 



(A.17) 


Aj(t) 

[ w dA / d(jj A ] 


YA 

2 c 2 


(A.18) 


where e A is the Alfven wave polarization vector. (Refer to page 36 for the 
relevant coordinate axes.) 


Equation (A. 15) thus becomes 
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Pa = 


lim _ 
T — ♦ oo 2 


J_ / d*k [ °du - 
n 2 T J ( 


s(u-u> A (k)) 41_V(A19) 

udA / dcj) h± 

' u=u>a 


lim 

T -» co 4 


/ **** j ° du s { u ~ ua ^) ~ • ^ ^- 2 °) 


fc, 2 


where w 0 is a cut-off frequency such that u 0 « 

Continuing in much that same way as in the calculation immediately 
preceding this one and using the Dobrowolny/Veltri expression for J( e \ we obtain 


Pa 


2 KVa ( K ‘ dk (dk, S mHk„L/2) S m^k,y o D/2V A ) . ,(k,D\ 
L k „ ‘ J kj {k<V„D/2V / ;)Hk 1 D/2j‘ V 2 / 


(A.21) 


where K 0 - u> 0 / V 0 . 

We can certainly assume ^«1. The integral over k y has been done 

in the preceding calculation. We can apply that result to obtain 


(A22) 


Since K 0 « 17 C ,- / V 0 , we can take K 0 D « 1 so long as D < 25 m. Then 
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sin(k x D/2) 

P/2) 


1 throughout the range of integration, which gives us 


r dkx (i n (^/2) 

A D C* J 0 V ' * x 2 


(A.23) 


Again we have arrived at the very same expression as before (equation (A. 11)) in 
the limit where k z D « 1 except for a factor of 2. This factor could be due to a 
mistake we’ve made, but the more likely source is from equation (A.16) which we 
inferred from the Dobrowolny/Veltri expression (A.l). 

These authors took an approach different from the one followed above. 
They attempted to carry through the complicated calculations using the general 
expression for A, Ay, and e, only applying the Alfven approximation at a later 
stage in the process. 

Having written down the general formula (A.l), the authors then proceed to 
give results based on calculations which they characterize as “lengthy,” too lengthy, 
evidently, to be summarized even in an appendix. Thus there is no way to follow 
their calculations step by step in order to judge this part of their work. They 
obtain (equations (13)-(15) in the original): 
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where 


P = ^ X / d * k j H {k,0,<p)G{k,O,<p,w)6(u - kM} 


dA 

du> 


.-i 


6[u - w CT (fc)] 


k^k^k^ 

n*z "ty "'z 


[A. 24) 


H(k,0,<p) = s\n 2 (k x ^sin 2 ^ |^sin 2 ^ 


(A. 25) 


G(k,0,<p,u>) = |^Cit 2 _ n^^cos 2 0 + eisin 2 0 ) 


. 2 [ c 3 + ft 2 sin 2 0]‘ 

+ 62 sin 2 ip - — 


cos 2 <p + 


c i c 2 - m 2 ^ 63 cos 2 0 + eisin 2 0 ^ 


(A.26) 


In formula (A.24), “the notation ^ denotes a summation over all plasma modes 

a 

which are possible in the cold-plasma approximation, u> = to c (k) being the 
dispersion relation for the a mode.” 

At this stage of their calculation, they insert the Alfven wave dispersion 
relation into their equations. 

There are two (presumably typographical) errors in the condition obtained 
for the angle 0 that k makes with the z axis. The authors quote 0 = 
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arc tan 


and say that V A / V 0 « 1 implies 6 « n/2. In fact, 6 = 


arc 


(- K -) 

\V A cos <pj 

tan( — — ), and the assumption (true for TSS) is V A /V 0 » 1. The 

\ y 0 cos <pj 


expression they obtain for the Alfven wave power is 


= 32 ^£ r kUk r . n0d0 i 

v A j 0 j 0 y_ r/ 2 u A {k) 

{H(kA<p) • G(k,0,<e,u A {k )) . [6(8 - tfO+ 6(8 - 0 2 )]} (A. 27) 


The Alfven dispersion relation is only valid for frequencies much less than the ion 
cyclotron frequency. Thus for the power calculation to make sense as a calculation 
of power into Alfven waves, where the Alfven dispersion relation is inserted into 
the integrand, the integral over u> in Equation (A.24) must be restricted to a range 
| u> | «n ci . We note that there seems to be another typographical error in 

equation (A.27). Presumably a factor [k x 2 k y 2 k z 2 dA/dcj] 1 should be included in 
the integrand. The more important point is that the integration over all w has 
been carried out and it has not placed any limitation on the size of k z in the 
integrand. If the rest of the integrand were such that only small values of k x 
were important, this would not be significant; but in the orbiting wire 
approximation, which is later carried further to an orbiting mathematical line 
approximation, this cannot be valid. 
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At this point the authors introduce the restriction on k z . The correctness 
of this step is more difficult to ascertain because of the likelihood of a 
typographical error in the preceding step. We did not attempt to reconstruct the 
calculations in full, since it would have required too many guesses about what 
expressions the authors were using for the dielectric tensor, etc. In any case the 
next step involves “very lengthy” calculations, of which no details are reported. 

It is only at this stage of their calculations that Dobrowolny and Veltri 
begin to make approximations about the system dimensions. Their approxima- 
tions are rather extreme and certainly do not correspond to a real tethered satellite 
system. Despite their statement near the beginning of the paper to the effect that 
the earlier work of Belcastro, Dobrowolny and Veltri [1982] had missed the Alfven 
wings entirely because it had assumed a tether of infinite length, they proceed to 
make the same approximation. Not only that— they also make the extreme 
approximation of an orbiting wire with negligible radius. This approximation is 
also at odds with their previous statement that the finite size of the system was of 
essential importance. The point is that, by making these approximations at this 
later stage in the analysis, the authors have made it difficult to judge their effects. 
The approximations do not seem to be consistent. At one point they take D = 
Dn ci /v 0 « 1 . Later they take ~ 1/D, which would imply £ m » 1, when the 
Alfven condition is « 1. 
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Using the D « 1 limit (orbiting wire) they had arrived at the result 
((35) in their paper) 


Pa = 



+ 



(A.28) 


They then take ~ -L . 

Now is the maximum value of k z V 0 /Q cl allowed for the Alfven 
dispersion relation to apply. It is thus necessarily much less than 1. Taking 

~ 4: violates that condition in the extreme. They then take £ m ~ 1 to obtain 
D 


Ra = 



(A.29) 


as the Alfven wave impedance. 

We note that they calculated the Alfven wave impedance at around lfi 
before this sleight of hand involving the use of £ m ~ 

Let us take their final expression (A.29) for this wave impedance and apply 
it to the 100 km long tethered system with a 1 mm diameter tether. 
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Since L/D = 10 8 , ~ 2 x 10 -4 , and corresponds to around 

0.060, we find R x =* 1200/7! 

Thus the authors seem to have made a serious mistake in obtaining their 
expression for the Alfven wave impedance. Perhaps they meant to apply ~ 
1/D only to the case for which D » 1 (though they should have said so). Then 
they would have obtained 


3 c 2 



L 

njD 3 


(A. 30) 


, 1 
(assuming — » 



But this expression is unlike any previous results and differs from that of 
i / yl 

Drell, et al. by a factor -I — ° 

3 \ V A lJ ci D 

Dobrowolny and less have not really argued that the quadratic dependence 
obtained by them is the correct one, only that a dependence on the velocity that 
clearly vanishes as the velocity goes to zero is necessary. There is no 

disagreement on this point. The time variation in the problem is strictly due to 
the relative motion between the satellite system and the magneto-plasma, so there 
would be no waves excited by a motionless tethered satellite, should such a system 
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be possible. As we shall demonstrate, the SAO calculation also gives a null result 
when the velocity is zero. It is also plausible that the Alfven wave impedance 
should not otherwise depend on the velocity, so that the “physical” argument for 
the Dobrowolny/Veltri result turns out to be weak. 

We have now arrived at the SAO expression three different ways, but 
let us examine the two parts of the criticism of the SAO results made by 
Dobrowolny and his colleagues. First, let us consider their assertion that the 
SAO results cannot be correct because they do not go to zero with the system’s 
velocity. This assertion is demonstrably false. Consider our expression for the 
Fourier transform of E y (equation (18) in our JGR paper). It contains the factors 
ujS(u - k x v x ). This is to be integrated over oj. Now the fundamental, defining 
property of the Dirac delta function is 



u> 0 )du = f(u 0 ) 


Thus when we take the inverse Fourier transform to obtain E y , which contains a 
linear factor of oj, we obtain zero if the satellite velocity is zero. All subsequent 
calculations thus assume that the satellite velocity is not zero, for indeed there are 
no waves generated in the zero velocity case. Thus for Dobrowolny and his 
colleagues to take us to task for failing to have an explicit dependence on the 
satellite velocity in all our subsequent results, which are various limiting cases (all 
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of which assume a non-zero satellite velocity), is unfair. 

less has maintained that, if one first does the inverse Fourier integration 
and then takes the limit as v x goes to zero, the SAO expression for the wave 
impedance does not go to zero. It is easy to demonstrate that one cannot expect 
always to get the same result by inverting the order of integration and limit 
taking. Consider 


f{u,x) = 


oj 6(<jj - u> 0 ) 

x 2 + u r } 


and take the integral over all u> and x. 


f f(u>,x)dxdu = f - - ^ — dx = tt, 

J J x 2 + <jJ 0 1 

which is independent of u 0 . But it has already been assumed that u> 0 is not zero, 
so that taking the limit as w 0 goes to zero after the integration is nonsensical. In 
physical calculations we must always be guided by the physics in deciding 
“mathematical” questions; as, for example, causality is frequently invoked to 
determine contour integration paths. Let us apply physical reasoning to the case 
at hand. What is the source of the u> factor in our integrand? It corresponds to 
a time derivative. What is the source of the S(u - k z v x ) factor? It comes from 
the physical condition that the only time variation in the problem is due to the 


motion of the system with respect to the plasma. Thus taking v x to be zero 
before integration is the choice that makes physical as well as mathematical sense. 

Finally, we note again that Dobrowolny and Veltri concentrate on explaining 
their differences with the results of Drell et al. presumably because both analyses 
found a similar dependence on the tethered systems’ dimensions. 

They do not discuss why the Drell formula’s linear variation with the tether 
length and inverse variation with the system’s dimension along the direction of 
motion should apply to the case they had considered at all, since Rasmussen et al. 
[1985] had already pointed out that the Drell results did not apply to a tethered 
system. In fact, Chu and Gross in 1966 ascertained that the Drell results were 
not applicable to any system having a vertical dimension that was not much 
smaller than its dimension along the direction of motion, i.e. just the opposite 
situation from what applies in the case of a tethered satellite. 
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